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We show that 1 // noise is produced in a 3D electron glass by charge fluctuations due to electrons hopping 
between isolated sites and a percolating network at low temperatures. The low frequency noise spectrum goes 
as co^" with a slightly larger than 1. This result together with the temperature dependence of a and the noise 
amplitude are in good agreement with the recent experiments. These results hold true both with a flat, nonin- 
teracting density of states and with a density of states that includes Coulomb interactions. In the latter case, the 
density of states has a Coulomb gap that fills in with increasing temperature. For a large Coulomb gap width, 
this density of states gives a dc conductivity with a hopping exponent of ~ 0.75 which has been observed in 
recent experiments. For a small Coulomb gap width, the hopping exponent ~ 0.5. 
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INTRODUCTION 

Low frequency 1 // noise |[ll0l2ll is found in a wide variety 
of conducting systems such as metals, semiconductors, tun- 
nel junctions |4jl, and even superconducting SQUIDs ||5l|6|]. 
Yet the microscopic mechanisms are still not well understood. 
One example is an electron glass which is an insulator where 
electrons are localized by a strong random potential. A spe- 
cial case of this is a Coulomb glass in which the electrons in- 
teract with one another via a long range Coulomb potential. 
Doped semiconductors and strongly disordered metals pro- 
vide examples of electron glasses. Experimental studies on 
doped silicon inversion layers have shown that low frequency 
1 // noise is produced by hopping conduction 1 7 1 . Because the 
systems are glassy, electron hopping can occur on very long 
time scales which can produce low frequency noise. In this 
paper we show that the resulting noise spectrum goes as 
where / is frequency and the temperature dependent exponent 
a > 1. 

Shklovskii has suggested that 1/f noise is caused by fluc- 
tuations in the number of electrons in an infinite percolating 
cluster 10] . These fluctuations are caused by the slow ex- 
change of electrons between the infinite conducting cluster 
and small isolated donor clusters. Subsequently Kogan and 
Shklovskii combined a more rigorous calculation with numer- 
ical simulations and found a noise spectrum where a was con- 
siderably lower than 1 (9|. Furthermore, below a minimum 
frequency of order 1-100 Hz, the noise spectral density satu- 
rated and became a constant independent of frequency. Their 
calculations were valid only in the high temperature regime 
where the impurity band was assumed to be occupied uni- 
formly and long-range Coulomb correlations were essentially 
neglected. Since then there have been attempts to include the 
effects of correlations. 

In particular, Kozub suggested a model lHoll in which elec- 
tron hops within isolated pairs of impurities produce fluctua- 
tions in the potential seen by other hopping electrons that con- 
tribute to the current. While leading to 1 //-type noise within 
some frequency range, this model also shows low frequency 
noise saturation due to the exponentially small probability of 



finding an isolated pair of sites with a long tunneling time. 
Moreover, the noise magnitude is predicted to increase as the 
temperature T ^ in contradiction with the recent experi- 
mental findings of Massey and Lee |11|. This, in part, led 
Massey and Lee to the conclusion that the single particle pic- 
ture is inconsistent with the observed noise behavior A dif- 
ferent approach was proposed by Kogan 1 12] who considered 
intervalley transitions as the source of the hopping conduc- 
tion noise. Unfortunately this approach does not seem to be 
analytically tractable and is not easily generalizable. 

In this paper we extend Kogan and Shklovskii's approach 
]9] by including the energy dependence of the hopping as well 
as the effects of electron-electron interactions on the single 
particle density of states ^(e). This is essentially a mean field 
approximation: we assume that charge is carried by electron- 
like quasiparticles whose interaction with the other charges 
is taken into account via the single particle density of states. 
Later we will present some justification for why we believe 
this approach works for low frequency noise. For comparison 
we also consider the case of noninteracting electrons with a 
flat density of states. 

The paper is organized as follows. In section Ila, we de- 
scribe our calculation of the noise spectral density. In section 
lib, we present the density of states that includes the Coulomb 
gap and that models the decrease in the gap with increasing 
temperature. We show that this form of the density of states 
yields the usual value of the hopping exponent 5 w 0.5 for 
small values of the Coulomb gap width E^. However, for large 
values of 5 ^ 0.75. Both values have been seen experi- 
mentally UH ElIIl 111 mill- In section III, we present 
our results. 



CALCULATION 

Noise Spectral Density 

We start with a model of the Coulomb glass in which elec- 
trons occupy half of the impurity sites. Each site can have 
at most 1 electron due to a large onsite repulsion. The sites 
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are randomly placed according to a uniform spatial distribu- 
tion, and each has a random onsite energy (|), chosen from a 
uniform distribution extending from — W/2 to W/2. Thus, go, 
the density of states without interactions, is flat. At T = such 
a system is a perfect insulator while at low but finite temper- 
atures it will be able to conduct via variable range hopping 
fTM fl?. 20]. In this picture the DC conductivity is dominated 
by particles hopping along the percolating network, which is 
constructed as follows. The resistance associated with a 
transition between sites / and j grows exponentially with both 
their separation r,y and energy difference 8,^: 



R,j=R°jexp{x,j) 



(1) 



where the prefactor Rfj = kT /{e^yfj) with yfj being given by 



,2|aJ| r2e2-|2 ^2 



3k^ 
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where D is the deformation potential, s is the speed of sound, 
p is the mass density, t, is the localization length and K is the 
dielectric constant. = Ej — e,- — e^/Kr,y is the change in 
energy that results from hopping from / to j with 8/ — (]),■ + 
Y,j 4p~nj being a single site energy. In Eq. ([1), the exponent is 
given by 



kT 



(3) 



The exponent reflects the thermally activated hopping rate be- 
tween / and /■ as well as the wavefunction overlap between the 
sites. 



max[|8,-~A'Ue;-A'|], - <") (E; ^ A*) > 
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(In what follows we choose the Fermi level /j = 0.) 

A noninteracting picture of DC conduction is described in 
terms of electron hopping between sites in a cluster that spans 
the entire sample. In order to determine which sites are in a 
cluster, we introduce the "acceptance" parameter x such that 
any two sites / and j are considered "connected" if Xij < x 
and disconnected otherwise. For small values of x only rare 
pairs of sites are connected. As we increase x, more such pairs 
appear and small clusters start coalescing into bigger ones un- 
til an infinite cluster - the critical percolating network - is 
formed at some Xc- At this point we can neglect the con- 
tribution of the remaining impurity sites to the DC conduc- 
tivity since it is exponentially small compared to that of the 
sites already in the percolating network (although the former 
sites are important for understanding both AC conductivity 
and noise). In the same spirit, the resistance of the critical 
percolating network is dominated by a few pairs with x,, — x^ 
- these are the pairs that bridge the gaps between large finite 
clusters enabling the formation of the infinite cluster. Hence, 



the resistance of the entire sample is well approximated by 
Rtot~R° exp (xc ) where ^ ee A:r / (e^yo) with y° being the av- 
erage value of Y°- given by Eq. 0. 

In the presence of Coulomb interactions, there is no exact 
mapping of transport onto a percolation picture. We never- 
theless assume that upon diagonalizing the interacting Hamil- 
tonian one finds that charge carrying excitations are of a lo- 
cal nature, and so they can be treated within the percolation 
picture as noninteracting quasiparticles. The Coulomb inter- 
actions renormalize the single-particle density of states which 
acquires a soft gap. We will discuss this in more detail in the 
section on the density of states. However, we will mention 
here that this approach appears to work well for DC conduc- 
tion and leads to a temperature dependence of the conductivity 
which is distinctly different from the noninteract- 
ing case and which agrees with experiment (see for example 
ref. I22I1). However, the question about the validity of this ap- 
proach is still far from being settled - see i24ll for a different 
point of view. 

In our treatment we will focus on the noise caused by quasi- 
particle hopping between isolated clusters and the percolating 
network, producing fluctuations of charge in the latter |8, ^. 
Let Np be the average number of such particles in the critical 
percolating network and dN^ (f ) be its time-dependent fluctu- 
ation. Assuming that only stationary processes are involved 
(i.e. {5N,j, {t2)?)N,j, (f 1 )) = f{t2 — 1\ )), we can use the Wiener- 
Khintchine theorem yl to relate the noise spectral density 
5/(co) of current fluctuations to the Fourier transform of the 
autocorrelation function: 



5/(co) _2(8jV^(f2)SjV^(fi))o 

/2 



Ni 



(5) 



where / is the average current. The charge fluctuation autocor- 
relation function can be expressed as a superposition of modes 
a, each of which relax exponentially with a characteristic time 
la- Thus the Fourier transform (. . .)(n of the autocorrelation 
function is a weighted sum over Lorentzians 01 • 



2kT 



(5Arj,(f2)5A^j.(fi))co = — ^ - 
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(6) 



Here C,- = {e^/kT^ fi{\ — fi) is the "capacitance" of site 
/ (with fi ~ [exp(8,/^r) + 1]^' being its equilibrium occu- 
pancy) while Tfj' and x|/a(0 are the a-th eigenvalue and eigen- 
vector of the following system of linear equations: 



(7) 



with Rij being the inter-site resistances given by Eq. ([Q. Since 
RJj^ is proportional to the hopping rate ijj^ — exp(— x,y) 
from site / to site j, eq. relates ijj^ to the relaxation rates 
Tfji' of the entire percolating network. The sum over sites / 
in Eq. (|6j runs only over those sites that belong to the critical 
percolating network (CN) since only their occupancies affect 
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the current through the sample. The physical meaning of the 
quantity C,i|/a(0 is that it is proportional to the fluctuation hfi 
of the occupation of site / and decays exponentially with the 
associated time constant T^. The eigenvectors satisfy the fol- 
lowing conditions: 



La c,¥a(0¥;(;) = 5u 



(8) 
(9) 

(10) 



The first condition states that the eigenfunctions are orthonor- 
mal; the second states that the functions form a complete set. 
One of the eigenfunctions is a constant which we take to be 
the one corresponding to a = 0. This has the eigenvalue 
Tq ' =0. Eq. dlOt is the orthonormalization condition be- 
tween this eigenstate and the others. It represents the fact that 
the fluctuations in occupation represented by the a 7^ modes 
do not affect the total number of electrons on the impuritiy 
sites. Thus the last equation is just the statement of overall 
charge conservation. We remark here that Eqs. are lin- 
ear only within the assumption made earlier of noninteracting 
quasiparticles. Otherwise the Rij are not constant coefficients; 
they depend on the onsite energies, which in turn depend on 
the occupancies of other sites. 

Since we are interested in the modes that affect the charge 
in the conducting network, we can replace the sum over a 
by a sum over all finite clusters that coalesce with the infi- 
nite cluster as the acceptance parameter increases above Xc- In 
particular we can replace the sum over a by an integral over 
X and a sum over all finite clusters merging with the infinite 
cluster at a given value of x. With this in mind, we can eval- 
uate Eq. (|6jl using Eqs. (|8} and dlOt . For a single mode a the 
sum over sites / can be split into a sum over finite clusters (FC) 
and a sum over the infinite cluster (IC). So we can write the 
normalization condition Eq. (|8|l and the charge conservation 
condition Eq. as 

C„ryifl{m)+Y^C„yifl{n) = \ (11) 

mGFC «GIC 

£ C^\|/a(m) + ^ C„^a{n) = (12) 
mGFC nGiC 

Since the fast modes equilibrate the occupations of sites 
within each cluster, the eigenfunctions do not depend on their 
site indices within each cluster, i.e., V(/a('«) = ¥a,FC, V m E 
FC and x|/a(n) — ¥a,ic, V m e IC. As a result the we can 
take i|/a.FC and \(/a.ic out of the sums in Eqs. ™d ( I12t . 
The sum over capacitances in the finite clusters will be much 
smaller than the sum over the infinite cluster which implies 
that (v|/a,ic)^ is negligible in Eq. (^Ji- This leads to 




-1/2 



¥a,FC 



Plugging this into Eq. il2\ yields 



¥a,ic Y Cm = - { C,ri 

melC VniGFC 
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(13) 



(14) 



We can use these results to evaluate the sum over sites in 
Eq. ^ by noting that all the sites in the critical network are 
also in the infinite cluster by definition. Thus 



1/2 



(15) 



where A^ic (x) is the number of sites in the infinite cluster at a 
given value of x. 

In evaluating Eq. (|6|l, we make the following approxima- 
tion for Xa- Since we are interested in the modes a that af- 
fect the charge of the percolating network, we only consider 
particle exchange between the isolated clusters and the in- 
finite cluster. This involves hopping times that are longer 
than those within the percolating network itself by defini- 
tion. Due to the exponentially wide distribution of hopping 
times Xij such exchange is likely to be dominated by the sin- 
gle closest pair of sites of which one belongs to the finite, and 
the other to the infinite cluster The relaxation times within 
each cluster are much faster, and therefore the above men- 
tioned pair serves as a "bottleneck" for intercluster relaxation. 
A simple diagonalization of the system of equations for 
two clusters Ai and J42, with the "bottleneck" hopping resis- 
tance R = mm{Rij; i E J4i, j E between them (and with 
the assumption that all other intercluster resistances are much 
higher and all intracluster resistances are much lower than R) 
leads to the foUowing expression for the intercluster relaxation 
time: 



x=R 



H -r 



EC, 

Ue^2 J 



-1 



(16) 



Since we are interested only in the situation when one of the 
clusters is infinite, this simplifies Eq. M6V . z ~ R J^ieji Ci, 
where is the finite cluster. 

We can substitute this value of T into Eq. (|6} by replacing 
the sum over all modes a by a sum over all finite clusters that 
coalesce with the infinite cluster as the acceptance parameters 
is increased above Xq. Each such finite cluster contributes one 
new term to the sum over a in Eq. with the corresponding 
la — T^iest Ci where R{x) — R°e^. Then we can write the 
spectral density of the noise as follows: 



16kT r N[c\x)R{x) (L'g^ CQ' 



(17) 



where stands for the sum over all finite clusters that coa- 
lesce with the infinite cluster as x increases by dx. The param- 
eter ^ > 1 and sets the distance in x space from the percolation 
threshold. 

This equation is difficult to evaluate mathematically. Fortu- 
nately, however, we can extract the low frequency asymptotic 
behavior of Eq. (I17> where the above approximations are well 
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justified. The lowest frequency contributions come from large 
values of x where the infinite cluster has already absorbed al- 
most all the sites (i.e. A^ic ~ A^, the total number of sites). 
What is left are the small clusters, which are mostly isolated 
sites in the increasingly rare voids of the infinite cluster. The 
probability of having two such sites in the same void is negli- 
gibly small. Since low frequency noise will be dominated by 
the hops between such isolated sites and the infinite cluster, 
we only consider such hops in obtaining the spectral density 
of current fluctuations. In Eq. we can set A, to correspond 
to this situation at large x, and we can replace the sum over all 
finite clusters that are merging with the infinite cluster with a 
sum over all sites multiplied by the probability Pi {x,t)dx that 
a single site with energy e has its nearest neighbor between x 
and X + dx. 

We can write down an expression for P\{x,t)dx. We be- 
gin by defining Pi(x,8) to be the probability that a given site 
with the onsite energy e has no neighbors nearer than x. Let 
exp [— p (-« , €)dx] be the probability that a site with energy £ has 
no neighbors between x and x + dx. Then 



Pi {x, e) — exp 



p(x',e)iix' 



(18) 



We can use this to express Pi (x, t)dx as the product of P\ (x, e), 
the probability of no neighbors within x, multiplied by the 
probability of having a neighbor between x and x + dx: 



Pi{x,E)dx = Pi(x,e) 



— Pi(x,e) 
ax 



dx 



(19) 
(20) 



Thus (— 3Pi/3x) is the probability density for a site to have its 
nearest neighbor between x and x + dx. We can now write the 
spectral density of current fluctuations as 



/2 



l6kTV r , f^'^ . , J dPi(x,E) 
- T / dx deg(e,T)[ -^-^ 

R{x)C^{€) 



(21) 



where V is the volume, W is the bandwidth, and /(e) is 
the Fermi occupation number. To obtain an expression for 
Pi (x, e), we note that the average number dN of impurity sites 
found in a phase volume element dQ. = d'^'rde! within a dis- 
tance jc of a site with energy e is given by 



dN = g{z')Q 



2r 



2kT 



^ ' ' de'd''r (22) 



The probability that no sites are in dQ. is given by 



energy 8 has no neighbors nearer than x is given by 

rW/2 



exp <— d r 



-W/2 



de'gie'j) 



2r |£| 
X B{x-j- — 



2kT 



(24) 



Notice the absence of the Coulomb energy in the argument of 
the 9-function in Eq. i24\ . in accordance with our quasipar- 
ticle picture. Our quasiparticle picture is likely to work best 
for hops between isolated sites and the infinite cluster. Al- 
though one such hop may result in a sequence of other hops, 
these will mostly happen within the infinite cluster on a much 
shorter time-scale, effectively renormalizing the properties of 
the "slow" particle. As was mentioned earlier, these renormal- 
izations can be included in the single particle density of states 

To facilitate evaluating the integral in Eq. J21l i numerically 
for the case where we include a Coulomb gap in the den- 
sity of states, we define the dimensionless variables r — r/t,, 
e = e/Eg,(h = &/f,T = kT/Eg, x = r^(jc)C(8) = /(e) (1 - 
/(e))e^, and |(e,f ) = g{e,T)/go. go is the noninteracting 
density of states and Eg « y^Kg^JT^ is the characteristic 
width of the Coulomb gap. Evaluating the integral over x in 
Eq. (I21> leads us to define 



(25) 



x~2r-\ 

Then we can rewrite Eq. ( I2U as 
5/(0)) _ r«'/2 



/2 



W/2 



d'egieX 



2T 



w/2 

de'g{e',T 

-w/2 ' Jo 



f^dr 



Q[x-KXc) 1 , ,-,9^9/,- (26) 



1 + a)2x2(x,e) 



where A ^ 6A%glElVi,^ / {N^'), Rv = (3y/47i)'/V^, W 
W/Eg, ri = AizgoEgi^, and 



Pi(i,e) 



exp 



6 



T'^dr' 



W/2 

I 

-W/2 



rfe"|(e",f) 



2f' 



|e| + |e"| + |e- 
2f 



(27) 



For comparison we also consider the case with no Coulomb 
gap by setting g{e,T) = go in Eqs. (I21> and i24\ . Since 
there is no natural energy scale, we do not rescale the en- 
ergies. However, we can define r, T, and cb as before. As 
a result, the definition of x in Eq. ( I25t becomes x = 2f + 
(|e| + |e'| + |e-e'|)/(2r). In Eq. (l26l. A is replaced bv 
Aa = 64nV g^^^ / N^y" and W is replaced by simply W. In 
Eq. i27\ r\ is replaced by r\o = 47t^^g„. 



^ dN 



-dN 



(23) 



Thus the probability Pi{x,e) that a given site with the onsite 



DENSITY OF STATES 

At zero temperature, long-range interactions produce a 
Coulomb gap centered at the Fermi energy in g(e, T) | IsiEll 
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l25l 12611 . This gap arises because the stability of the ground 
state with respect to single electron hopping from an occu- 
pied site ; to an unoccupied site j requires that the energy 
difference A;* > 0. At finite temperatures the Coulomb gap 
is partially filled and the density of states no longer vanishes 

at the Fermi energy 113 S El EH El ■ The exact form 

1 I 



of g{£-,T) is not known, but some have argued 130. l3lL l3^ 
that its low temperature asymptotic behavior is described by 
^(e = 0,r) ~ r''^^ We have done Monte Carlo simulations 
of a three dimensional Coulomb glass with off-diagonal dis- 
order and we find that ;?(e — 0,T) cannot be described by a 
simple power law |28, 33]. The results of such simulations do 
not produce a density of states that is suitable for use in our 
noise integrals due to finite size effects. In particular g{e, T) 
goes to zero at energies far away from the Fermi energy be- 
cause of the finite size of the system. 

Another way to approximate the density of states is to use 
the Bethe-Peierls-Weiss (BPW) approximation |31|. The 
idea is to treat the interactions between one "central" site and 
all other sites (boundary sites) exactly, but to include the in- 
teractions between these boundary sites by means of effective 
fields. The density of states can then be written as a convolu- 
tion 



Wo/2 1 / p' 

de'g(e~e') y-hi^ 



(28) 



where ^(e) is the zero temperature density of states and Wg is 
the bandwidth. The function h [z/kT) takes into account ther- 
mal fluctuations in the occupation of the central site and the 
boundary sites. At low temperatures it has a sharp peak with 
a width of the order kT s.iz = 0. We can make the approxima- 
tion (1 /kT)h (e/kT) w -/'(e) where /'(e) is the derivative of 
the Fermi function. The zero temperature density of states can 
be determined numericaUy by solving a self-consistent equa- 
tion based on the ground state stability condition that a single 
electron hopping from an occupied site ; to an unoccupied site 
j requires a/ > lE3.Eil- The result of evaluating Eq. ( I28> 
is shown in Fig.[3 

Since using the BPW approximation to evaluate Eqs. i2l\ 
and J24t is rather awkward, we model the finite temperature 
density of states by 



{kTf 



{kT) 



2 • 



(29) 



Notice that for T = 0, g{e, T = 0) ~ e^ for e < Eg as is ex- 
pected for a Coulomb gap in three dimensions. For large ener- 
gies (e ^ Eg and e ^ kT), g{e,T) approaches the noninteract- 
ing value go. A comparison of Eq. (I29> with the BPW approx- 
imation at various temperatures is shown in Fig.^ Eq. ( I29I I is 
the expression we use for the density of states of a Coulomb 
glass in Eqs. (I21> and i24\ . 

We can calculate the DC conductivity resulting from this 
density of states by following Mott's argument for variable 
range hopping 1 18]. We start with the hopping resistance Rij 
given by Eq. Q. Mott pointed out that hopping conduction 
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FIG. 1 : The density of states g(e, T) versus e at various temperatures. 
The symbols are calculated using Eq. with Wo/2 = 2.3 x 10"* K. 
The density of states is measured from the Fermi energy Ef = 0. 
The lines are the result of evaluating Eq. j29j with Eg — 100 K. go — 



6.25 X 10 states, 



s/KA^ 



at low temperatures comes from states near the Fermi energy. 
If we consider states within e„ of the Fermi energy (Ef = 0), 
then the concentration of states in this band is 



N{eo,T)= I gie,T)de 



(30) 



where g{e, T) is given by Eq. i29l . So the typical separation 
between sites is /?„ ~ [N{£o-T]^^^^ . To estimate the resistance 
corresponding to hopping between two typical states in the 
band, we replace with /?„ and e,; with eo in Eq. (|3j to obtain 
x{eo). Minimizing x(eo) numerically yields e,^. Aplotofx(eo) 
versus temperature is shown in Figure |2| The dc conductivity 
is then given by (y{T) = aoexp[— x(eo)]. We find that at low 
temperatures (T <^ Eg) 



o{T) = OoCXp 



(31) 



where 5 is the hopping exponent. The value of 5 depends on 
Eg. For large values of the Coulomb gap (Eg 



'50K)5; 

< 



iO.75 

while for small values of the Coulomb gap (Eg ^ 1 K) 5 w 0.5. 
When we tried intermediate values of T — 8, 10, and 20 K, 
we found that ln[x(eo)] versus ln(r) had a break in slope with 
5 ~ 0.5 at low temperatures and with 5 w 0.72 — 0.75 at high 
temperatures. Examples are shown in Fig. |5] 5 = 0.75 is 
higher than the Mott value of 5 = 0.25 associated with a flat 
density of states and the value of 5 = 0.5 derived by Efros 
and Shklovskii |21] for the zero temperature Coulomb gap. 
However, experiments on materials such as ultrathin metal 
films find values for 5 = 0.75 ± 0.05 El H [H E El 
in agreement with our value of 5 for large Eg. The mecha- 
nism behind this exponent has been a puzzle Here 
we see that a possible simple explanation for the experimen- 
tal observation of an anomalous hopping exponent is that the 
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FIG. 2: x{Eo,T) versus temperature with go = 6.25 x 10^^ states/K- 
and ^ = lOflo = 5.29177 A where cio is the Bohr radius and 10 is 
an estimate of the dielectric constant. We show plots for Eg = 0.4 K 
(O), Eg = SK (□) and Eg = 200 K (A). Eg = 0.4 K corresponds to 
the value of the Coulomb gap deduced from transport measurements 
while £■(, = 8 K value from tunneling measurements on Si:B 1 22, 37]. 
The lines are fits to the numerical data with the indicated slopes. The 
fit to the Eg = 0.4 K data yields 5 = 0.56 and = 19 K. The fit to the 
£g = 8 K data at low temperatures yields 5 = 0.47 and Tg = 21, 206 
K, while the fit to the high temperature data yields 6 = 0.72 and 
5 = 0.357 K. The fit to the Eg = 200 K data yields 5 = 0.75 and To = 
42, 068 K. 5 is virtually independent of go but To does depend 

go ■ 

For example, changing go by 10 orders of magnitude to 6.25 x 10+^ 
states/K-i^ results in 8 = 0.75 and To = 19 K for Eg = 200 K. 




FIG. 3: The noise power spectrum as a function of frequency. The 
frequency is measured in the units of Y° which is estimated to be of 
the order lO'^Hz for values appropriate for insulating Si:B. Unless 
otherwise noted, all curves in this and the following figures which 
were obtained for the case with a Coulomb gap used T| — A-TlEgbj gQ 
12[£g/(e^/K^)]^ = 4.8 X 10^'', which in our estimates corresponds 
to the experimental dopant concentration of roughly n = 0.8;Jc for 
Si:B 111, 23]. We set W = 20, Ry = 100 and hcc = 1 (the precise 
value of X has no effect on the low frequency noise which is governed 
by X ^ Xc). The parameter A = 64lzVEgg^t,^ /N^y°. For comparison 
we show the noise spectrum in the absence of a Coulomb gap with 
g{£,T) = go in Eqs. <21> and <24> . In the absence of a Coulomb 
gap, A is replaced by A,, = 647tV g^t^^ / N^y" and r| is replaced by 
^o = 47l^Vo = 4.8 X lO^*'. The energy is measured in arbitrary units 
and we set W = 20. The other variables are the same as in the case 
of a finite Coulomb gap. 



Coulomb gap in the single particle density of states is fill- 
ing in with increasing temperature. If one takes this into 
account in the variable range hopping calculations, then the 
observed exponent of 0.75 can be obtained naturally. How- 
ever, we should caution that our calculation applies to three 
dimensions while a two dimensional calculation may be more 
appropriate for ultrathin films. In fact we find that the anal- 
ogous two dimensional calculation with a density of states 
gie,T)^ go m+kT)/ (Eg + |e| + ^r) yields 5 « 0.5. 

RESULTS 

We evaluate Eqs. i26\ and i27\ numerically and display the 
results in Figs. |3}]6l In Fig.|3]we show the spectral density 
of the noise as a function of frequency. We find that for a 
wide range of parameters the noise spectral density is given 
by S{03) ~ CO^" with the spectral exponent a between 1.07 
and 1.16 (see Figs.|3]|5} which is "1/f" noise. For comparison 
we show in Fig. |3lthe noise spectrum in the absence of a 
Coulomb gap with §(e, T) = go in Eqs. (12 1> and (I24> . The 
slope of a line through the open squares is — 1 . 12 which is very 
close to the values obtained with a Coulomb gap. Notice that 
the presence of a Coulomb gap reduces the noise amplitude at 
low temperatures. 



In Fig.|3]we use the transport value of Eg sa OAK, not the 
tunneling one ~ S/T; the two were found to be different by an 
order of magnitude l2^ E3ll . We find that increasing Eg by a 
factor of 20 does not produce a noticable change of the results 
at low temperatures (T = 0.1 Eg), but at high temperatures 
(T = 10 Eg) it does lead to saturation of the noise power at 
low frequencies. This is shown in Figurel^which also shows 
that saturation occurs in the absence of a Coulomb gap when 
Tjo is increased by a factor of 20. This saturation of the noise 
power occurs because the probability Pi {x, e) of finding a site 
with no neighbors closer than x (see Eq. ( I24» decreases expo- 
nentially with increasing temperature and with increasing r| or 
r|„. In addition Pi {x,e) becomes exponentially small as x be- 
comes large, and it is the large values of x that contribute to the 
low frequency noise. Finally we note that decreasing Eg by a 
factor of 10 does not produce a noticable change of the results 
for either low temperatures (T = 0. 1 Eg) or high temperatures 
(r = 10 Eg). We plot the spectral exponent a in Fig. Inversus 
temperature for the cases with and without a Coulomb gap in 
the density of states. In both cases we see that it decreases 
slightly with increasing temperature and eventually saturates 
in qualitatively agreement with experiment ■ Fig.|6lshows 
that the noise amplitude \/S grows with temperature and even- 
tually saturates, both in good qualitative agreement with the 
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FIG. 4: The noise power spectrum as a function of frequency at T = 
10 for various values of r| = AllE^i^go. The rest of the parameters 
are the same as in Fig. |3| Notice the saturation at low frequencies for 
large r|. For comparison we show the case with no Coulomb gap at 
r = 10 with a large value of r|o = Ani^go. Large values of r\o lead 
to saturation but small values do not. 



experimental results of Massey and Lee fTHl . The data of 
Massey and Lee span 2 decades in frequency while our cal- 
culations are able to cover a much broader range. Again we 
see from Fig.|6lthat the presence of a Coulomb gap reduces the 
noise amplitude at low temperatures. We obtain qualitatively 
the same results both with and without a Coulomb gap in the 
density of states which implies that the behavior of the noise 
spectral density with respect to temperature and frequency is 
not strongly tied to the hopping exponent 5 or to the particular 
form of the density of states. 

We will now discuss some of the physical reasons behind 
our results. The fact that we obtain 1// noise is perhaps to be 
expected since weighted sums over Lorentizians (see Eq. (|6j) 
often result in 1 // noise 1 1]. The subtlety lies in the tempera- 
ture dependence of the noise amplitude. For simplicity let us 
consider the case of a density of states with no Coulomb gap 
which gives qualitatively the same results as the case with a 
Coulomb gap. The decrease in the noise amplitude \/S with 
decreasing temperature is due to the presence of activated hop- 
ping processes which decrease with decreasing temperature. 
However, this is not at all obvious from Eq. M6\ . The integral 
for the noise power at low frequencies is dominated by large 
X which corresponds to long relaxation times T ~ exp(i). In 
this case the factor of /(e) [1 — /(e)] cancels between the nu- 
merator and denominator leaving the temperature dependence 
of the integrand dominated by Pi(x,e)exp(— i). Pi(x,e) in- 
creases while exp(— i) decreases with decreasing temperature. 
The fact that our calculations yield an decrease in the noise 
amplitude with decreasing temperature implies that the acti- 
vated hopping processes associated with exp(— i) dominate. 
We should mention that experimentally the noise power does 
not always decrease with decreasing temperature. In some 
cases it increases with decreasing temperature L38>..391 but we 



FIG. 5: The spectral exponent a as a function of temperature with 
a Coulomb gap in the density of states (o) and with a flat density 
of states (O)- We have suppressed the error bars for the case with 
no Coulomb gap to avoid cluttering the graph. The suppressed error 
bars are comparable to those for the exponent with a Coulomb gap 
at high temperatures. The temperature is measured in units of the 
Coulomb gap Eg for the case where there is a Coulomb gap, and in 
arbitrary units for the case without a Coulomb gap. The inset shows 
the experimental data obtained for Si:B | 111 . 



10' 
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2 4 6 8 10 

Temperature 

FIG. 6: Noise amplitude \/S at O) = 10" '^y" (or / ^ IHz) as a func- 
tion of temperature for the cases with a Coulomb gap (o) and without 
a Coulomb gap (Q)- The temperature is measured in units of Eg for 
the case of a finite Coulomb gap and in arbitrary units in the case of 
no Coulomb gap. The inset shows the experimental data for / = IHz 
liill. 



do not know the differences in the samples which can account 
for this difference in behavior. 

To summarize, recent experiments on 1 // noise flUl are 
consistent with a quasiparticle percolation picture of transport 
in electron glasses, though this does not exclude multi-particle 
correlations. 
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